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ABSTRACT 


The  reported  research  presents  results  on  the  fracture  strength 
of  graphite/epoxy  (T300/5208)  [0/±45]  and  [0/+45/90]  symmetric  lami- 
nates in  the  presence  of  angle  cracks.  The  results  indicate  that  the 
fracture  strength  can  be  predicted  irrespective  of  the  crack  orienta- 
tion if  one  uses  the  effective  normal  crack  length  in  the  equation 
developed  for  normal  cracks.  The  fracture  strength  prediction  is 
made  using  a two  parameter  (unnotched  fracture  strength  and  inferred 
dimension  of  the  damage  zone)  model,  when  the  laminate  is  subjected 
to  uniaxial  loading  along  one  of  the  material  symmetry  axes.  The  re- 
sults indicate  that  the  critical  damage  zone  is  not  constant  for  the 
two  laminates  although  it  is  independent  of  the  crack  orientation. 
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INTRODUCTION 

A review  of  the  literature  indicates  that  most  of  the  work  in  fracture 
of  composites  has  emphasized  Mode  I loading.  By  comparison,  mixed-mode 
fracture  has  received  little  attention.  Most  of  the  available  mixed-mode 
data  falls  in  the  category  of  unidirectional  composites,  with  the  crack  being 
parallel  to  the  fibers.  For  this  case  it  has  been  observed  that  cracks 
oriented  parallel  to  the  fibers  propagate  along  the  fiber  direction  [1]. 

In  addition,  Wu  [1]  showed  that  the  critical  stress  intensity  factors  ^Ic 
and  are  practically  independent  of  crack  length  for  fiber  glass  rein- 
forced Scotchply. 

For  co-linear  crack  propagation  a constant  maximum  energy  release 
rate  was  suggested  [1]  as  a possible  fracture  criterion.  However,  the 
data  for  Scotchply  indicates  that  this  criterion  is  not  realistic. 

McKinney  [2],  using  unidirectional  graphite/epoxy  laminates  with 
implanted  cracks  parallel  to  fibers,  also  found  that  the  energy  criterion 
was  not  applicable.  Sanford  and  Stonesifer  [3]  reported  contradicting 
results  for  glass  reinforced  plastics  which  seemed  to  show  the  validity 
of  the  energy  criterion  although  the  data  does  not  include  a measurement 

'^IIc* 

All  of  the  mixed-mode  fracture  data  previously  cited  was  generated 
for  unidirectional  composites  with  crack  parallel  to  the  fiber  direction. 

By  comparison,  less  is  known  about  composites  where  the  layup  is  more 
general  than  the  unidirectional.  Brinson  and  Yeow  [4]  studied  fracture  of 
graphite/epoxy  laminates  using  specimens  cut  at  various  angles  to  the 
principal  material  direction.  Their  results  indicate  that,  in  general,  net 
fracture  stresses  tend  to  increase  when  the  crack  is  aligned  more  closely 
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with  the  load  direction.  No  attempt  was  made  to  analytically  model  the 
fracture  behavior  of  single  and  double  edge  notched  specimens. 

Some  of  the  available  criteria  for  mixed-mode  fracture  of  isotropic 
materials  were  discussed  by  Nuismer  and  Hahn  [5].  On  the  whole,  all  the 
criteria  except  the  one  based  on  the  energy  release  rate  assume  that  the 
final  fracture  becomes  imminent  when  failure  occurs  at  a point  some  dis- 
tance away  from  the  crack  tip.  The  failure  functions  considered  include 
the  circumferential  stress  [6],  the  principal  stress,  the  Tresca-Guest  and 
Huber-Mises-Hencky  yield  functions  and  a modified  version  [7],  and  the 
strain  energy  density  [8].  Recently,  Wu  [9,  10]  proposed  the  use  of  the 
polynomial  failure  criterion  to  predict  fracture  of  glass/epoxy  laminates. 

In  the  present  work,  however,  no  attempt  will  be  made  to  apply  any 
of  the  aforementioned  methods  to  the  fracture  of  graphite/epoxy  laminates. 
Rather,  results  are  presented  on  the  variation  of  the  notched  strength 
with  the  crack  orientation  angle  relative  to  the  load  axis.  A simple 
model  based  on  the  idea  of  an  effective  normal  crack  will  be  chosen  to 
fit  the  data.  This  model  can  be  useful  in  assessing  the  severity  of  bal- 
listic impact  damage  in  composite  laminates. 


EXPERIMENTAL  PROGRAM 


The  experimental  program  consisted  of  determining  the  unnotched 
and  notched  fracture  strengths  of  two  composite  laminate  systems.  It 
will  be  shown  that  a two  parameter  model  can  be  used  to  predict  the 
fracture  strength  of  the  laminates  when  subjected  to  uniaxial  loading 
along  one  of  the  nnaterial  symnetry  axes. 

The  material  used  in  the  experimental  program  consisted  of  Thornel 


300  graphite  fibers  in  a Narmco  5208  epoxy  resin.  All  test  specimens 
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(see  Figure  1)  were  2-in.  wide  and  12-in.  long  (9-in.  between  end  tabs). 
Cracks  were  produced  by  first  drilling  a small  hole  in  a specimen  followed 
by  a final  lengthening  with  a 5 mil  diamond  wire.  No  attempt  was  made  to 
further  sharpen  the  crack  tips.  Two  crack  lengths  were  used,  0.5-in.  and 
0.8-in.  The  angle  0 (Figure  1)  varied  from  15°  to  90°,  in  15°  increments. 

The  laminate  orientations  and  number  of  test  specimens  may  be 
Summarized  as: 

[0/±45/90]g;  three  specimens  for  each  crack  length  and  angle 

[0/±45]g;  three  specimens  for  each  crack  length  and  angle 
Thus,  a total  of  72  specimens  were  tested  to  determine  fracture  loads 
as  a function  of  laminate  orientation,  crack  length,  and  angle  of  crack. 

The  results  are  given  in  Table  1.  In  addition,  four  specimens  of  each 
laminate  orientation  were  tested  to  determine  the  unnotched  tensile 
strength;  the  results  may  be  found  in  Table  2.  The  fiber  volume  fraction 
for  each  laminate  was  approximately  0.64. 

All  specimens,  both  notched  and  unnotched,  were  loaded  by  friction 
grips  and  tested  in  a closed-loop  hydraulic  machine  at  a constant  load 
rate  of  20  Ib/sec.  The  environmental  conditions  were  approximately  con- 
stant at  70°F  and  50%  relative  humidity. 

RESULTS 

Several  models  have  been  proposed  in  the  literature  for  the  prediction 
of  the  fracture  strength  of  composite  materials.  Two  models  that  are  of 
interest  in  regard  to  the  data  herein  are  the  inherent  flaw  model  [11]  and 
the  average  stress  model  [12,  13].  Both  models  have  been  used  to  predict 
the  fracture  strength  of  uniaxially  loaded  composites  containing  holes 
and  cracks.  In  the  case  of  cracks  the  average  stress  model  [12,  13]  becomes 
identical  with  the  inherent  flaw  model  of  [11],  as  discussed  in  [14]. 
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According  to  the  above  models,  the  ratio  of  the  notched  fracture 
strength  to  the  unnotched  tensile  strength  is  given  by 


f ^0  \ 
K I 


1/2 


(1) 


where: 


noo 


'0 


notched  fracture  strength  for  an  infinite  plate 

unnotched  tensile  strength 

original  crack  half-length 

dimension  of  the  damage  zone  at  crack  tips 


The  model  given  by  Egn.  (1)  has  been  applied  to  composite  laminates 
with  center  cracks  perpendicular  to  the  load  axis.  It  will  be  shown  that 
the  same  model  can  be  applied  to  the  angle  crack  problem,  where  the  angle 
crack  is  treated  as  an  equivalent  normal  crack  having  a length  equal  to 
the  projection  of  the  angle  crack  normal  to  the  load  axis.  In  other 
words,  the  original  crack  half-length,  ag,  is  replaced  by  the  equivalent 
normal  crack  half-length  of  magnitude  ag  sin  0,  where  0 is  the  angle 
between  load  axis  and  crack  (Fig.  1).  The  damage  zone  dimension,  Cg, 
then  retains  the  same  meaning  as  for  a normal  crack. 

The  infinite  plate  notched  fracture  strength  is  calculated  from  the 
equation 


a = Ya 
n“  n 


(2) 


where:  Y = 


finite  width  correction  factor  based  on  equivalent 
normal  crack 


a_- 


notched  fracture  stength  of  finite  width  specimen, 
based  on  gross  cross  sectional  area 


The  value  of  Y is  determined  from  the  equation  for  isotropic  materials  [15], 
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Y = 1 + 0.128  (a/w)  - 0.288  (a/w)^  + 1.525  (a/w)^  (3) 

where,  for  the  angle  crack  problem,  2a  is  the  equivalent  normal  crack 
(2a^sin  0)  and  2w  is  the  specimen  width.  The  finite  width  correction  fac- 
tor for  the  [0/±45]g  laminate  is  slightly  different  from  the  isotropic 
value.  However,  the  difference  is  less  than  Z%  in  the  range  of  crack 
sizes  considered  [16]. 

The  value  of  the  damage  zone,  Cq,  is  found  empirically.  Rewriting 
Eqn.  (1)  in  the  form 

results  in  an  equation  of  a straight  line,  with  Cq  as  the  slope.  Note  that 
Eqn.  (4)  has  been  written  for  the  angle  crack. 

Figure  2 shows  the  results  of  the  [0/±45/90]^  laminate  with  angle 
cracks.  Three  specimens  of  each  crack  length  and  angle  were  tested. 

Using  a least  squares  fit  the  value  of  Cq  was  found  to  be  0.113  in.  (2.87  mm). 

It  is  of  interest  to  compare  the  value  of  Cq  found  using  angle  cracks 
with  that  found  using  only  normal  cracks,  for  the  same  material  system. 

The  normal  crack  results  are  shown  in  Fig.  3,  for  the  oata  of  Ref.  [17]. 

The  values  of  Cq  for  angle  and  normal  cracks  are  within  6%  from  the  average 
value. 

Figure  4 compares  the  model  prediction,  Eqn.  (1),  for  both  the  angle 
crack  and  normal  crack.  Note  the  insignificant  difference  in  results  re- 
gardless of  the  type  of  crack.  Also,  one  additional  comment  is  in  order. 

The  data  for  the  normal  crack  is  for  a [0/90/±45]g  laminate  [17].  However, 
the  data  of  Whitney  and  Kim  [18]  indicates  that  the  notched  strength  of 
quasi -isotropic  graphite/epoxy  laminates  appears  to  be  independent  of 
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[ 

; stacking  sequence.  Thus,  for  the  stacking  sequence  indicated  in  Fig.  4, 

I use  of  either  the  normal  or  angle  crack  results  in  essentially  the  same 

\ value  of  predicted  notched  strength.  A comparison  of  Eqn.  (1)  with  ex- 

i 

I perimental  results  will  be  shown  in  a later  figure. 

! To  be  of  use,  the  value  of  Cq  in  Eqn.  (1)  must  remain  constant  for 

Iall  crack  sizes.  Figures  5 and  6 show  the  effect  of  crack  size  on  the 
damage  zone  dimension,  for  angle  and  normal  cracks,  respectively.  The 
^ solid  lines  represent  the  least  squares  value  of  Cq.  A constant  value  of 

ICq  for  different  size  cracks  appears  adequate. 

Corresponding  results  for  [0/±45]  laminates  may  be  seen  in  Figs.  7- 

f 

[ 11.  All  normal  crack  data  was  taken  from  Ref.  [17].  Figure  7 shows  data 

[ for  angle  cracks,  whereas  Fig.  8 combined  data  for  two  laminate  thicknesses. 

[ The  effect  of  thickness  within  the  framework  of  the  classical  laminated 

i 

; plate  theory  is  illustrated  in  Fig.  11,  which  shows  different  values  of 

i 

I Cg  for  the  two  laminates.  Since  smaller  values  of  Cg  means  higher  notch 

' sensitivity,  i.e.,  more  reduction  in  notched  strength,  the  [0/±45]g  laminate 

r 

is  seen  to  be  slightly  more  sensitive  to  cracks  than  is  the  [0/±45]25  lami- 
nate. However,  the  difference  is  practically  negligible  and  the  average 
value  is  used  in  the  subsequent  analyses. 

Figure  9 compares  best  fit  fracture  strengths  for  angle  and  normal 
cracks.  The  normal  crack  curve  is  drawn  using  Cg  from  Fig.  8.  As  before, 
it  is  seen  that  predicted  results  are  essentially  the  same  regardless  of 
whether  Cg  is  found  using  angle  crack  data  or  normal  crack  data. 


Figure  12  shows  both  experimental  results  and  analytical  predictions, 
Eqn.  (1),  for  two  different  laminates  with  angle  cracks.  The  results  indi- 
cate that  the  analytical  model  adequately  predicts  the  notched  strength  of 
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laminates  under  uniaxial  loading  conditions.  Even  though  there  is  more 
scatter  in  data  for  the  quasi -isotropic  laminate,  the  analytical  model 
gives  results  that  are  conservative  on  the  average. 

For  the  [0/±45]  laminates  the  value  of  Cq  obtained  using  the  least 
square  fit  method  is  less  than  the  value  of  Cq  obtained  using  compliance 
matching  [17].  Consequently,  the  compliance  matching  method  resulted  in 
strength  prediction  higher  than  the  experimental  data. 

Table  3 summarizes  the  average  values  of  Cq  for  the  laminate  systems 
and  crack  types  studied.  The  results  of  this  table  indicate  that  the 
critical  damage  zone  size  Cq  is  not  constant  for  different  laminates  although 
it  is  independent  of  the  crack  orientation.  Thus,  Cq  should  be  found  for 
each  laminate,  as  was  done  in  the  present  study. 

CONCLUSIONS 

It  has  been  shown  that  a two  parameter  (unnotched  fracture  strength 
and  inferred  dimension  of  the  damage  zone)  model  can  be  successfully  used 
to  predict  the  fracture  strength  of  composite  laminates  subjected  to 
uniaxial  loading  along  one  of  the  material  symmetry  axes.  The  model  is 
applicable  regardless  of  the  orientation  of  the  crack,  or  the  laminate 
system  (at  least  for  those  laminate  systems  that  have  been  tested).  For 
the  material  system,  T300/5208,  and  laminates,  [0/±45]g  and  [0/±45]2g, 
the  data  are  essentially  independent  oF  laminate  thickness. 

However,  the  value  of  the  dimension  of  the  damage  zone  is  different 
for  each  laminate  system.  Table  3 gives  a comparison  of  the  values  of  Cq 
for  the  laminate  systems  and  crack  types  studied.  The  results  of  this 
table  indicate  that  the  damage  zone  Cq  is  not  constant  for  various 
laminates  although  it  is  independent  of  the  crack  orientation. 
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Table  2.  Ultimate  Strengths 


Laminate 

Width 

Thickness 

Strength 

(in) 

(in) 

(ksi) 

[0/±45]g 

0.9462 

0.0351 

71.66 

0.9920 

0.0351 

70.65 

0.9477 

0.0344 

72.70 

0.9864 

0.0335 

64.46 

Average 

0.9681 

0.0345 

69.87 

C.V.  (%) 

2.39 

4.92 

5.21 

[0/±45/90]g 

0.9495 

0.0453 

58.12 

0.9908 

0.0457 

59.63 

0.9497 

0.0427 

55.51 

0.9908 

0.0436 

54.40 

Average 

0.9702 

0.0443 

56.91 

C.V.  (%) 

2.45 

5.03 

4.47 

’ 
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Table  3.  Critical  Damage  Zone  Sizes 


Laminate 

Type  of  Crack 

Cq,  in.  (mm) 

[0/±45]5 

Normal 

0.054  (1.372) 

[0/±45]2s 

Normal 

0.054  (1.372) 

[0/90/±45]g 

Normal 

0.100  (2.540) 

[0/±45]5 

Angle 

0.049  (1.245) 

[0/145/90]^ 

Angle 

0.113  (2.870) 

M/CaujsOo) 
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Figure  3.  Least  squares  fit  for  normal  cracks, 
[0/90/±45]5  laminate.  ^ 
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Figure  7.  Least  squares  fit  for  Cq,  angle  crack,  [0/±45] 
laminate. 


Figure  10.  Damage  zone  dimension  vs.  equivalent  normal  crack  for  [0/+45]  laminate. 
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Figure  11 


Effect  of  crack  length  and  thickness  on  Cf,  for  [0/+45] 
laminate,  normal  crack.  ^ 


